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I. INTRODUCTION 



Catastrophes of all kinds can be roughly defined as short duration-large amplitude events fol- 
lowing and followed by long periods of "ripening". Major earthquakes surely belong to the class 
of 'catastrophic' events. Because of the space-time scales involved, an experimental approach is 
often difficult, not to say impossible, however desirable it could be. Described in this article is a 
"laboratory" setup that yields data of a type that is amenable to theoretical methods of prediction. 
'^C^" ' Observations are made of a critical slowing down in the noisy signal of a solder wire creeping under 

I ' constant stress. This effect is shown to be a fair signal of the forthcoming catastrophe in both of 

, two dynamical models. The first is an "abstract" model in which a time dependent quantity drifts 

. slowly but makes quick jumps from time to time. The second is a realistic physical model for the 

collective motion of dislocations (the Ananthakrishna set of equations for creep). Hope thus exists 
that similar changes in the response to noise could forewarn catastrophes in other situations, where 
. such precursor effects should manifest early enough. 

^ , Catastrophes as defined in the abstract are related to a class of phenomena sometimes called "relaxation" oscilla- 
te ■ tions. The observation of two widely separated time scales makes relaxation oscillations a good a priori subject of 
^ I theoretical investigation, because one may suspect that their formation results from the existence of a (more or less 
^ ' hidden) small parameter. This gives some hope of a "general" theory based upon the small size of the parameter. 
. ^ To take an example, some major earthquakes lasting a few tens of seconds occur in the same general area about 
^ . every hundred years. They thus involve a ratio of 'typical' times in the neighborhood of 10~^, which is a very small 
number. In the present study a laboratory experiment was devised which displays similar relaxation oscillations, but 
acting time-wise on a scale that is convenient for investigation. Experimental observations are explained in light of 
a two-part dynamical bifurcation model of catastrophe. Comprising the local form of a physical model shown to be 
valid when close to the catastrophe, the striking result is as follows. In response to an external source of noise the 
two models (the local one and the physical one) predict fluctuations with a correlation time that increases before the 
catastrophe, just as is observed experimentally. The important point is that this well known critical slowing down 
phenomenon occurs significantly before the transition and could be used to forewarn it. 
If^ We study the creeping of a soft metal under constant stress. Accurate time records show that this creeping actually 

displays the following time dependent noisy component. The wire typically lengthens slowly with background (non- 
thermal) noise, with sometimes a "large" sliding event, followed again by a noisy slow lengthening regime, etc. 

Plastic deformation of solids is a complex phenomenon, not yet fully understood [l|. It has long been observed to 
take place in a non-smooth manner. Most studies have focused on the Portevin-Le Chatelier effect observed under 
constant strain rate conditions, whereas the present experiment is concerned with creep at constant stress. In both 
cases large steps (cyclic slips) are embedded in a noisy background. There is a general agreement that this is due to 
the complex dynamics of networks of dislocations whose motion is a means for the stressed solid to flow. As reported 
in section [n] below, careful observations of the creep in strained Sn-Pb quasi-eutectic material (a solder wire) show 
the following : 

i) On average a sample under constant stress lengthens at constant rate. 

ii) Continuous monitoring shows time dependent fluctuations of this length superposed on its secular increase. 
Hi) From time to time the length jumps by steps. Afterwards a noisy and steady (on average) length increase is 

recovered until the next jump, etc. 

Explaining all this remains a challenge for the common models of creep. Nevertheless it is of great interest because 
it can be seen as a laboratory model of other far less accessible phenomena like earthquakes, where on average there 
is also a continuous slow sliding with random microseismic noise, interrupted by large fast sliding steps characterizing 
major earthquakes. We recently introduced the idea that in such systems the slow to fast transition can be 
described by a saddle-node bifurcation in a dynamical system evolving slowly with time. We pointed out the interest 
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in such a modelization, which allows one to theoretically predict the response of this dynamical system to a noise 
source. It was shown that the induced fluctuations drift toward low frequencies (i.e., toward large correlation times) 
before the transition, and so could be used as a forewarning. 

Here we show that the dynamical saddle-node bifurcation model is the reduced form of a set of equations previously 
used to describe the Portevin-Le Chatelier effect in metals, or metallic alloys, i.e., the Ananthakrishna (AK) model 
The saddle-node model describes fast transitions resulting from the intrinsic dynamics of the original system 
which can be described locally (close to the step) as a slowly rocking potential system. Using the AK model as well 
as the saddle-node reduced model we show that, with an added external source of noise, the correlation time of the 
fluctuations increases before the transition, following the classical scenario of slowing down at bifurcation points. As 
was shown long ago by Dorodnytsin the same local dynamics describe the transition from slow manifold to fast 
transients in relaxation oscillations of dynamical models, like the van der Pol equation in the strongly nonlinear limit. 
By looking at the fluctuations of the length of our samples, we found this behavior near the step-like transitions, with 
the characteristic drift to low frequencies before the transition. 

Unlike recent publications that have presented the idea [3| that precursors of earthquakes could be found in the 
mechanical response to external perturbations (such as the increase of the fluctuations and their slowing down near 
transitions) our idea goes further, by giving an order of magnitude of the precursor time. The authors of do not 
introduce the effect of a given time dependence of the parameters and consider only systems with steady parameter 
values on both sides of the bifurcation. To quote reference Q "The suggested approach to analytical study of any 
kind of catastrophes is based essentially on the solution of a stationary problem of the possibility and conditions of the 
unstable equilibrium state in the system in question" . Without taking into account explicitly the time dependence of 
the parameters sweeping the bifurcation set, it is impossible to get the time scale for predictions. As we show, this 
scale depends crucially on the rate of change of the parameters near the bifurcation, which may be estimated from 
the knowledge of the ratio of the two time scales (fast and slow ones) for saddle-node models. 

An important challenge lies in the difhculty of stating a suitable model for a given catastrophic event. Actually 
there is more than one class of possible slow-to-fast transitions in dynamical systems. The dynamical saddle-node 
can be seen as belonging to the class of systems with an equilibrium point losing stability as a parameter changes. 
It is not a loss of stability, but rather a loss of existence of the equilibrium point, occurring at the folding point of 
the slow manifold. However slow-to-fast transition may happen without any folding of the slow manifold. As shown 
in exploring a very often used mathematical representation of stick-slip behavior, the Dieterich-Ruina equations, 
the slow-fast transition can originate from the finite time singularity of the slow dynamics itself. In that case, critical 
speed-up, or drift of fluctuations toward large frequencies, is found to replace the critical slowing-down effect. This 
nonlinear phenomenon is obviously outside the class of phenomena explainable by a stability analysis of equilibria 
of dynamical systems. We refer the interested reader to the paper on this subject @. In the present paper we do 
not consider this case, because it is clearly not the one observed in the creeping experiments. Note that at this time 
it is unknown if real earthquakes (as well as other observed catastrophes) belong to the saddle-node case with a 
slowing-down near the transition or to the finite time singularity case with speed-up expected. 

The creeping experiment is described in section |TT] together with the striking spectral observations. In the two 
sections following section [Hi we present our theoretical approach. We introduce in section Hill an abstract dynamical 
model showing fast jumps. In this model, one assumes that, as a dynamical system, the jump follows a "saddle- 
node" bifurcation. There a pair of fixed points, one locally stable, the other locally unstable, merge and disappear 
as a control parameter changes. Afterwards, the system moves quickly to a new equilibrium state that is at finite 
distance (in phase space) from previous equilibria, whence the jump. After reviewing the saddle-node bifurcation 
in this light, we assume that the parameter changes with time, namely that the parameter defining the bifurcation 
is itself a slow function of time. When, by this change, the parameter crosses the bifurcation value, the dynamical 
system makes an abrupt transition and jumps "generically" from one equilibrium state to another. In section HVl we 
show that this "abstract" model is pertinent for describing the slow-fast transition in the relaxation regime of a set 
of equations derived by Ananthakrishna for creeping in the relaxation regime. This set of equations describes the 
dynamics of populations of dislocations in the creeping solid. The model shows relaxation oscillations in which slow 
drift is interrupted by fast variations. Near the slow-fast transition, we show that this model reduces to the generic 
equation quoted above, for a certain range of parameters. Therefore it displays a typical critical slowing-down in its 
response to external noise, in agreement with the real data reported in section |TT1 Of particular interest is the fact 
that, from the experimental data, one can predict in advance a "large slip" event. The event is preceded by a shift 
toward low frequencies in the random fluctuations of specimen length. 
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II. CREEPING EXPERIMENT 



Choose a uniform "wire" that creeps under small stress at room temperature, fabricated from a much studied 
material, Sn-Pb solder, alloyed to be nearly eutectic. To establish creep at virtually constant stress, let one end of 
the wire be fixed and at the other end establish a constant force of tension. By this means we found the length of 
the wire to increase on average at constant rate. By highly accurate monitoring of this length vs time we observed 
a small time dependent part with the following pattern: on a background of fluctuations, from time to time a large 
slip event was observed, after which the continuous lengthening with a small background noise was recovered. The 
observed noise is non-thermal, since thermal noise has a far too small amplitude to be relevant. It is assumed 
to result from rearrangement of defect structures in the poly-crystalline structure of the quasi-eutectic, involving 
dislocation dynamics. We assume that this noise originates from a source that is to first-order independent of the 
overall lengthening of the wire. This is equivalent to saying that it is due to ongoing micro-scale events triggered 
by the imposed stress, independent of the global creeping. Therefore we analyzed the response to this noise source 
according to the AK equations subject to an external noise (see below). 



A. apparatus 



The instrument used in these experiments, which is pictured in Fig|l] is an extensometer jSj]. Young's modulus 
can be accurately measured with a wire specimen, by placing different size masses on the weight pan. The trace of 
the vertical wire holding up the boom/weight pan in Fig [T] has been enforced (colored in black) to be visible in the 
image. Though not presently used, the black clamp was for purpose of holding a power resistor that was employed to 
measure the specimen's thermal coefficient of expansion. To measure temperature changes, a solid state thermometer 
was placed down into the sample space, through the hole seen near the top knurled clamp. 




FIG. 1: Extensometer used in the experiment. 



To calibrate the instrument a tungsten wire of diameter 0.1 mm was mounted in the extensometer, and signal 
output level changes were recorded as various gram-mass standards were placed on the weight pan. By using the 
known Young's modulus for tungsten, the resulting measurements yielded a constant of 1.0 nm per analog to digital 



count, for the 24 bit adc employed, which is sold by Symmetric Research [11|. This constant is applicable to the data 



presently reported. A different measurement technique yielded essentially the same calibration constant. A He-Ne 
laser was used with a mirror, operating as an optical lever, to measure boom position change as different masses were 
placed on the pan. 

On the basis of two factors, ordinary solder was chosen for the present study. First of all, an unusual property of 
tin is well known, when an ingot of the metal is strained by large amounts. The sound which it then emits is well 
described by the German word "zinngeschrei" , which translates "tin cries" . 

The initiative for this work was also influenced by the observation of unusual spectral features in the output from 
a novel seismo grap h [l2l | . The unusual low-frequency motions of the Earth's crust that were then observed to precede 
an earthquake 'l3l| are readily seen by the VolksMeter. This is due to the instrument's use of a "displacement" sensor, 
rather than the "velocity" sensor used by conventional seismometers. It was therefore natural to consider an alloy of 
tin, with the expectation that its defect properties should be more like those of the earth than is possible for a pure 
metal. The tin alloy for our study was ordinary soft solder (60% Sn/ 40% Pb), used universally in the electronics 
industry. Although the pure 'Eutectic' alloy is actually 63% Sn/37% Pb, we will nevertheless use this word to describe 
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our specimen in the discussions that follow. 

The stress level due to the load placed on the Pb-Sn wire used in the present experiments was considerably smaller 
than the typically lOMPa used in usual creep studies. The present load was due solely to the weight of aluminum 
comprising the boom plus empty weight pan of the extensometer. This stress value was estimated to be O.bMPa, 
based on considerations that include the distance of the wire's attachment point from the position of the (fine) roller 
bearing in the unpright housing, that supports the boom on its end opposite the pan. 



raw' creep data 




FIG. 2: (a) Elongation of the wire versus time, raw signal x{n). (b) Filtered signal y{n) at the exit of the single pole high-pass 
filter p. 

During a typical avalanche, the elongation of a 13.5cm long, 1.6mm diameter specimen, would be about 30/im. 
Before the avalanche, typical lengthening velocity is 5fim/s and typical rms fluctuations observed in the wire length 
(with secular term removed by high-pass filtering) would be about 50nm. 

An example of creep record (length growth versus time) is shown in Figl21[a). The fluctuations of the raw signal 
are not visible because they are much smaller than the average length variations. 



B. Data analysis 

To eliminate the average growth, we use a standard technique of filtering. At the exit of a low-pass single pole 
filter, the filtered signal y{n) is given by the recursive formula 



yn= ^ {Xn - Xn-l) + PVn-l, (1) 

where x{n) stands for the raw creep signal. In other words the filtered signal is the convolution product of the 
raw signal derivative by an exponential response function R{t) — exp— 27r/ct, where fc is the corner frequency of 
the filter. In equation ([Ij the parameter p is given hy p — exp— 27r/c(5t where St is the sampling time of the record, 
equal to j^sec in the experiment. Using a corner frequency fc — 50mHz, the filtered signal takes the form shown in 
Fig|2jb). In order to analyze the spectral properties of the fluctuations, and see how they evolve with time, we must 
consider separate sequences of a given time duration. One significant difficulty in this analysis is the proper choice 
for the magnitude of this time duration. The window cannot be too large, lest the spectrum becomes invariant with 
time. Conversely, if the window is too narrow, there will be excessive noise in the spectral (or correlation) signal. We 
believe the time windows have been judiciously chosen; so that there is a significant increase of the correlation time 
before a "burst" . Such increase of the correlation time (actually the decrease of the frequency width of the signal) 
was found to happen before every sliding event, giving some hope that this "critical slowing down" is pertinent for 
predicting the "catastrophe" before it occurs. For the present work a time duration of order 1/10 th of the interval 
between adjacent bursts was chosen. The spectral density of a given sequence (n^, Ui + N) of the filtered signal 



, ni+N 

^(") = ^ 2/(fc)exp2i7rnfc|2, (2) 

changes significantly during the creep process. The striking effect is the shift of the spectral density toward low 
frequencies in the last stage of the slow regime, i.e., just before the burst. Examples of this phenomenon are presented 
in the figures ([3])(a-b), which show the spectra corresponding to the time intervals (1,2) indicated in Figure ©-a. 
The two spectra clearly differ. The first spectral density (a) displays a large number of components; whereas the 
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FIG. 3: (a) and (b) Normalized spectral density of y{n) before the second step (the step starts ai n — 26370), for the two 
time intervals (1-2) indicated in Figure ©-a, at the abscissa 18370 < n < 20370 and 24370 < n = 26370 respectively, each 
for a time interval of A'' = 2000 counts, (c-d) Cumulative spectra corresponding to figures (a-b) respectively (e)Width w{t) of 
the cumulative spectrum versus time (at height equal to 75 per cent of the maximum), the vertical arrows indicate the very 
beginning of the fast steps. 



second spectral density, corresponding to the time interval just before the burst (figure b), is concentrated close to 
zero frequency. To quantify this effect, we calculate the cumulative spectrum 

Cs(n) = X]^W. (3) 
1 

which is a smooth curve whose asymptotic value (for n=N) gives the experimental variance of the fluctuations during 
the sequence considered 

= Cs{N). (4) 




FIG. 4: (a) Creep data versus time, (b) Width of the cumulative spectra before the burst with TV — 2000. The horizontal 
segments ending at each point indicates the time interval (before this point) over which the width is calculated. 

The shift of the density spectrum toward low frequencies is visually more clear in its associated cumulative spectrum, 
as seen in figures (c-d). We define the "low frequency extension" w of a sequence by the width of the cumulative 
spectrum corresponding to 75 per cent of its maximum value. 
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Cs{w) = ^CsiN). (5) 

Shown in figure 3(e) is the evolution of the spectral width for the whole experimental record. A drop in w is clearly 
seen to occur before each burst. 

The important point is that the decrease of the width before the burst is experimentally foreseeable because it 
occurs during a "precursor" time which is noticeably larger than the duration of a sequence N6t. This is illustrated 
in FigHJb) which displays a zoom of the spectral width evolution before the burst. 

The two following sections are devoted to theoretical models, the first one shows the main features of what is 
measured in the experiment on creeping, namely the critical slowing down of the fluctuation spectrum occurring 
before the sliding event, the second one is a model of creeping found in the literature which consists in a set of 
nonlinear coupled ordinary differential equations, derived by Ananthakrishna and collaborators [Tsj . 



III. DYNAMICAL SADDLE-NODE BIFURCATION 



This section explains why a saddle-node bifurcation with a slow sweeping of the bifurcation parameter exhibits a 
slowing down in its response to a source of noise in a window of time extending well before the bifurcation itself. 
This "abstract' model has no direct connection with the physical phenomenon of creeping. In the section afterwards, 
however, we explain that a model of creeping shows the same slowing down in a range of parameters, linked to a local 
saddle-node bifurcation. 

Consider first the saddle-node bifurcation of a "gradient flow", that is a damped dynamical system such that a 
coordinate x{t) is a solution of the equation of motion of the form 



dx dV 
dt = 

In this equation V{x) is a potential, and the dynamics tends to everywhere lower the value of V{x). 

The catastrophe theory of Thom and Arnol'd [lo| studies how steady equilibria of equations like ^ change under 
smooth deformations of the potential V{x). Below we consider a different kind of question, namely what happens to 
the solution of equation ([6]) when the potential V{x) becomes itself a slowly varying function of time, and particularly 
when a pair of equilibrium points disappears by a saddle-node bifurcation. Indeed this question of the sweeping across 
bifurcations has been already widely studied [l4| with various applications in mind (by " sweeping" we mean crossing 
of a transition point with a time dependent parameter in the equation(s) of motion). However, to the best of our 
knowledge the occurrence of an intermediate time scale in the case of slow sweeping has been overlooked, although 
we believe it to be crucial for a strategy of foretelling catastrophes in the real world. 



A. local cubic potential 

The equation ^ is too general to be very helpful. However, it may describe a saddle-node bifurcation where a stable 
equilibrium disappears, assuming that V depends slowly on time in a prescribed way, to become a function V{x,t). 
Near the transition, one may use a mathematical picture which is correct for a short time around the transition if the 
potential V{x,t) is a smooth function (see below for what happens afterwards). 

Assume first that V{x) does not depend explicitly on time and takes the form 

Vix)^-ilx' + bx), (7) 

with b real constant (for the moment). 

For b negative V{x) has two real extrema (i.e. the roots of ^ = 0), one —y/—b is a stable equilibrium, the other, 
y/—b, is an unstable equilibrium. For 6 = the two equilibria merge and disappear for b positive, see FiglSja). 
This is the saddle- node bifurcation. The shape of V{x) near x — and for b small is universal: for a given smooth 
V{x) showing this saddle- node bifurcation, one can always rescale x and the external parameter to obtain the "local' 
problem in this form. 

The extension to a time dependent control parameter b goes as follows. If & is a smooth function of time, one can 
assume that b{t) crosses the critical value, i.e. zero in the present case, at time zero in such a way that b{t) = at+ ... 
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with a non zero constant, dots being for higher terms in the Taylor expansion of h{t). For t and x close to zero, 
after rescaling, one can represent the dynamical system ([B]), close to the saddle- node bifurcation, by the "universal" 
parameterless equation 



Ax 9 



(8) 




(a) -3t \ (b) 

FIG. 5: (a) Cubic potential for h — —1,0, 1. (b) Quartic potential, b = —1,0, 1 



Outside of the neighborhood of x = 0, the solution of (jH) depends on other parameters defining V{.) for finite values 
of explained below. 

Although this is not obvious from its formulation, this model is also valid for the transition from slow to fast motion 
in the van der Pol oscillator in the limit of large non linearities, as shown by Dorodnitsyn 0|. Let us sketch the proof 
of this (interesting) point. Actually we shall look at the formally more general situation of relaxation oscillations, 
namely at solutions of a set of coupled ODE's (Ordinary differential equations) with a large parameter in the form: 



and 



X = r]F{x,y), 



V = G{x,y). 



(9) 



(10) 



In this set of equations dots are for time derivatives and the functions F{) and GO are smooth with values of order 
1 when their argument is also of order one. Moreover, is a large parameter. The slow manifold is defined by the 
condition that, in the limit rj large, the function F(x, y) must be close to zero over at least part of the trajectory. The 
Cartesian equation F(x, y) — defines a curve in the plane {x, y) which allows to find y as a function of a;, at least 
locally. This defines the equation of motion (along the slow manifold) 

y = G{x{y),y), 

where x{y) is such that F{x{y),y) = 0. The slow trajectory so defined stops at "folds" where the function x{y) ceases 



to be well defined, namely for values of {x, y) such that ^ 

discrete set of points. Near those points, the equation of motion can be solved by taking 5x 
where (a;o,2/o) are the Cartesian coordinates of the point such that 



— F^y = 0, F{x, y) = 0, F^yy being not zero. This defines a 

xo and Sy = y-yo 



F^F^,^ 0. (11) 

Let us look at the solution of the coupled equations ([1(11) near (xq, yo). The equation (ITUl) is not singular at this 
point and so can be solved for small variation of y and for t small, like 6y — tGa, t — being the (arbitrary) time 
where the trajectory is in O and G{xo^yo) — Gq. Consider now the first equation ([5]), and expand its right-hand side 
for 6x and 6y small. Because F^x = at {xq, yo), the first nontrivial term in the Taylor expansion of i^(.) near {xo,yo) 
with respect to x is ^{F^xx)oSx'^ . On the other hand the first term coming from the expansion with respect to Sy is 

F^ySy = (Fy)oGoi. Therefore, for t and 6x small, the equation for 6x derived from ([9]) reads: 

Sx^Tj Q (F,,)ofe' + {F^y)oGot] , (12) 
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One can check that all terms not written explicitly there are effectively negligible compared to the ones kept. Standard 
rescalings allows to transform the equation (jl2p into the " universal" equation ([5]) , provided various constraints of sign 
are satisfied by the quantities F^^x, F,y and G{.) all computed at {xo,yo)- This is a short version of the classical 
calculation by Dorodnycsin, showing that the " universal" equation ([5]) is also relevant for the slow to fast transitions 
in relaxation oscillations. It is worth mentioning also that this derivation does not really assume that the overall 
motion is periodic, as it may be extended quite easily to a system of many more coupled ODE's with only one fast 
variable. This shows that such large jumps are also possible in non periodic dynamics, because the dynamics on the 
slow manifold may be chaotic in between the jumps if this slow manifold has a number of dimensions sufficiently large. 

Now we shall focus on the slowing down near the dynamical saddle- node bifurcation described by equation (jH) . We 
shall first give its explicit solution. We look for a solution of equation ([5]) transiting from the "stable" fixed point at 
"large" negative times to the rolling down towards positive value of x at positive times. This solution behaves like 
x{t) « —\J~^ at large negative times. The equation is of the Riccati type and can be integrated by introducing 
the function such that = — ^ where y = ^ and y{t) is a solution of Airy's equation y + ty = 0. 

The solution of equation ([S]) relevant with the given condition at i — )■ — oo is drawn on Figure (|B]-a). In terms of the 
variable y{t) it is the Airy function Ai{—t) which writes 



Y{t) = Ai{-t) = / cos(— - ut)du. 
Jo 3 

Yet we have only solved the transient problem near the saddle-node bifurcation. The transition ends-up when t 
becomes equal to the first zero of the Airy function Ai{—t), i.e. the smallest root of the equation 

Y{t) = 0, a pure number, about tc ~ 2.338. It corresponds to a divergence of x{t) — — which behaves as 



just before this transition, as derived by expanding Y{t) close to tc. 

Therefore the "generic" equation ^ for the dynamical saddle- node bifurcation displays a finite time singularity. 
Let us precise the following mathematical subtlety. This property of the local flow, which results from the folding 
of the slow manifold, differs qualitatively from the finite time singularity found in the Dieterich-Ruina equations @ , 
where it was a property of the flow reduced to the slow manifold which is everywhere convex, as discussed in the 
introduction. In the case of the dynamical saddle-node bifurcation the singularity requires one to consider both the 
dynamics on and off the slow manifold, and this happens because the geometry forbids the continuation of an exact 
trajectory on this folded slow manifold. Actually the solution (fT3|) loses its physical meaning sometime before the 
singularity since the "universal" dynamical equation was derived under the assumption that x remains close to 
zero. This local theory cannot deal with finite variations away from the critical values, therefore we shall need to add 
finite amplitude effects to limit the growth of the instability after the transition, see the subsection IIII Bl 

We shall study now two questions, first the response of this dynamical system to an external noise, then the 
dynamics of a system showing a saddle-node bifurcation of the type just studied and reaching a new stable fixed point 
after this bifurcation. We explore first the response of our system to a small external noise, and look for qualitative 
changes in this response which could be a signal that occurs before the transition. 

Let us consider the equation ([8]) with a small noise added, so that equation ^ is replaced by 



x^x^ +t + e^(t), (14) 

where ^(t) is a random function of time, and e a small coefficient. 

In the limit e small, one can solve equation (jl4[) by expansion in powers of e, x{t) — XQ{t) + exi(t) + .... 
where 



The linear response to the noise is 

1 



Mt) = yj^ I dtat) Y^t). (15) 
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Because tends rapidly to zero as (<) tends to minus infinity, one can take <o = — oo to get rid of the effect of 

the initial conditions. 

To make the developments above more concrete, let us take a delta-correlated (or white) noise, such that < 

^{taMh) >=S{ta-h). 

The pair correlation of xi{t) is given by 

I rinl{t,t') 

where inf (t, t') is the smallest of the two real numbers t and t' . The behavior of this pair correlation for large negative 
values of both t and t' , is derived from the asymptotic expression of Airy's function, Ai{—t) ~ 2^[~'ty/^ ' Setting 
w — and F{w) 1 — w'^/'^ , the variance of the fluctuations writes 
< x,{tf >« (-i)/r ^iHef(-*)^^^^^W. 

In the limit (— i) — > oo the integral is concentrated near w = 1 so that 



<x^{tf>^-{-t)-"^, (16) 

which shows that the amplitude of the fluctuations increases some time before the transition itself. As the transition 
approaches, the standard deviation of the fluctuations, (j{t) = < {x{t) — xo{t))^ >, increases close to the critical 
time tc, because Y(tc) = 0. 



B. Quartic potential 



As the zeroth order solution diverges at < = tc, it does not make sense to describe the dynamical behavior of the 
fluctuations due to the external noise very close to as shown above. 

As said before, this unbounded growth of the fluctuations is a consequence of the local cubic form of V{x) when 
expanded near a; = 0, in obvious contradiction with the fact that x{t) tends to infinity. To suppress the divergence 
of x{t) after the saddle- node bifurcation we add a stabilizing (positive) term to the potential V{x) which becomes 
quartic, 



V,{x)^-^-bx+^, (17) 

as drawn in FiglSj^b). Because of the growth of Vq{x) at infinity, like x^, the solution of the differential equation 

x = b + x^-x^. (18) 

does not diverge at finite time. Notice that, formally this does not apply to cases where the local equation ([8]) describes 
a slow-to-fast transition in a limit cycle, instead of a gradient flow dynamics. However, unless one insists to look at 
what happens well after the saddle-node bifurcation, the details of the dynamics after the fast slide are not signiflcant, 
because one can consider that the point of landing on the slow manifold after the fast drift is like a new equilibrium 
for a gradient flow, this neglecting the slow motion on this manifold. 

The equation (jl8p can be written in the given scaled form for any quartic potential provided the coefficient of x^ is 
positive. For such a potential one parameter only remains. In equation (jl7p . we choosed to keep as explicit coefficient 
the coefficient b of the linear term in equation (fT7|) . For 6 = the dynamical system (jl8p is exactly at the saddle- node 
bifurcation, because at 6 = x = both the first and second derivative of Vq (x) vanish, but not the third derivative. 
Contrary to the case of the pure cubic potential, this system has always, that is for any value of 6, a stable fixed point 
beyond the pair of fixed points collapsing at the saddle-node bifurcation. This makes it a fair candidate for describing 
the dynamical saddle-node bifurcation without blow-up. 

As in the previous case, we shall take now a time dependent &, that will be taken as b — at with a positive constant. 
Because of the rescaling of the cubic and quartic term, the parameter a cannot be eliminated (another possibility 

3 4 

would be to put a parameter in front of the cubic term). For the potential Vq{x) — — o,tx -I- ^ we shall analyze 
the solution of the dynamical equation 

„2 ^3 



x = at + x-x-^, (19) 
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tending at large negative and positive times to the equilibrium point x = {aty^^, t being considered as a parameter, 
see FigElJa). Moreover we consider the limit a small, a being related to the ratio of small to large time scales, it can 
be estimated from experimental data. We shall prove that in this limit a small, there are three characteristic time 
intervals, depending how x is close to zero. 

The long time scale is the time lapse between successive major slips, typically of order 150 — 300s in our experiments 
on creeping. In our model it is the time needed for the potential Vq{x, t) to change significantly, to move from a pair 
of fixed points to a saddle-node bifurcation. Because time enters in Vq{x,t) through the combination (at), the a- 
dimensional time needed for a change of shape of Vq is of order 

icreep ^ ^ . (^0) 

The short time tsUp is of order unity in units of our model equation (jl9p as shown in the next paragraph. It is the 
duration of the abrupt change of slope of the function of time x{t), the rising of x{t) at tc- In the creep experiment 
reported above, it is about Is. Therefore the ratio of these two time scales tsUp/tcreep is as small as 10"^ — 10^^ in 
our experiments. 

There is another time scale, to, the time interval standing before the transition, and close to it, during which the 
potential is very flat, while the solution has not yet jumped. During this time, x and at are much smaller than unity, 
then the cubic term on the right-hand side of equation (jl9l) is negligible. In this range one recovers the universal 
equation of the dynamical saddle- node bifurcation ^ by taking X = xa~'^^'^ and T — to^/"^, with the boundary 
condition X{t) fa —T at T tending to minus infinity. This property concerns the rectangular domain drawn on 
Fig. |Bl-(a), where x is small, x ^ a}!'^ , and t extends from —a'^l'^ Xo t ^ a~^/^, located before the abrupt increase. 
Therefore the time extension of this domain introduces the intermediate time scale, 

to ^ (21) 

long compared to the short time t^up (which is of order unity, see below) and small compared to t creep = oT^ the 
average time between slips. 

We now show that the short time is of order unity, by matching the solution X(T) of the universal equation to the 
solution of equation ([2^ below, in the vicinity of the critical point tc(a) = a~^l^tc- Because X(T) behaves like ^ 
before it diverges, it follows that the solution x(t) behaves as « a-^^^t -t ^'^^ "large" values of St = a^^^'^tc — t before 
the critical time. This becomes of order one when St becomes also of order one. When this happens, the term at in 
equation (|19p is negligible, therefore the solution of this equation, which can be matched with the solution near the 
bifurcation, is the solution of the integrable equation 

i^x'^-x^, (22) 

with the asymptotic behavior for very large negative times x{St) ~ — This equation shows that, in our model, the 
fast time scale is of order one, because it has no explicit dependence with respect to the small parameter a. This 
result is confirmed by the numerics, see FigEJa). More precisely defining the rising time of x{t) by 1/5 of the width 
of the time derivative i, we get 

tslip ^ I7 (23) 

independently of the value of a, for a small. In the experiment, the intermediate time scale is tp'*^* ^ ^"^^^'^^up" ' 
which gives about 5s for a = 10~^, t^'^^^" = Is. 

In summary , by matching the two solutions in the range 1 ^ (—St) <^ a^^^"^, we show that the catastrophe takes 
place during this time tsUp which is of order one, because the displacement is then of order one, compared to the 
small displacement of order a^^^ taking place during time to = a~^^^ typical of the "universal" transition process. 

From this understanding of the various scales in the deterministic part of the dynamical equations, we can now 
look at the response to noise of this system, particularly at the range of time where something like a " critical slowing 
down" could be observed, and which actually happens in our experiments. 

With a noise source added, the dynamical equation © becomes, 

X = x^ - x^ + at + e^{t). (24) 

Actually the effective noise amplitude is not equal to e close to the saddle-node, as it depends on the value of the 
parameter a. Indeed for |t| < to , the cubic term in equation (I^Sj) is negligible, and the equation reduces to 

X = x'^ +at + eC{t). (25) 
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Xo 




FIG. 6: (a) solutions of equation (|14|) . with and without noise, e = (smooth curve) and large amplitude noise e — 1. (b) 
Solution of equation (|19|l for a = 10~^. The rectangle around the origin defines the region —to < t < to and —1/to < x < 1/to, 
with to ^ a~^^^ . The critical time is tc ~ 2.34 to- The two vertical lines inserted between the two arrows delimitate the large 
slope time duration, of order unity, (c) Solution of equation (p5)l with large amplitude noise e = 1 



which may be written as ^ = X"^ + T + e{a)^{t), by setting X = xa"^/^, T = ta^/^, and i{a) = ea"^/^. Therefore 
the effective noise is larger than e, by a factor a^^^^ in the rectangular domain of figEl 

Let us consider now the fluctuations of the solution x{t) of equation ((25)) . For a small noise source, the solution 
may be expanded in power of e as above. At first order it gives 

ii = [2xQ{t) - 3xl{t)]xi{t) + C{t), (26) 

whose solution is formally 

xi(i)= / diai)exp[g{t)-g{i)]. (27) 

Jto 

In general g{t) is the time integral of the second derivative of the potential — ^"i^^ , which yields with our choice of 

v{.y. 

git) = / [2xa{u) - 3xl{u)]. 



J to 

The standard deviation Gx-^ (t) has to be calculated numerically. We expect it to display the same behavior as for the 
cubic case in the whole domain where x{t) << 1 , i.e. a little before the transition and close to it, because the potential 
is cubic in this range. After the transition, we expect that the fluctuation decreases, because the solution without 
noise becomes quasi-steady. This is confirmed by the numerics; the amplitude of the fluctuations strongly increases 
close to the critical time tc, its maximum occurring at time tc, then it decreases. More precisely the standard deviation 
behaves exactly as xo{t), the red curve in Fig.(I3a), for small noise. Therefore the strong increase of the variance of the 
signal fluctuations cannot be used as a precursor for predicting the transition because it occurs simultaneously with 
the signal itself close to the transition. Note that this observation seems to contradict the currently found statement 
that fluctuation enhancement precedes the transition and can be used as a precursor. In the case of the saddle-node 
bifurcation model, we have indeed observed a "precursor" growth of the fluctuations, but only in the case of "large" 
amplitude noise, see FigjSt^b). Let us focus on the case of small noise. In this case we show below that the correlation 
time of the fluctuations changes much earlier than the onset of amplitude growth. This slowing-down can be used to 
foretell the event itself, in considerable advance of its occurrence. It is seen to consistently happen in the creeping 
experiment described below. 

Consider the case of small effective noise, where the correlation function and the spectrum of the fluctuations 
(x(t) — XQ{t)) are well described by the correlation function and spectrum of xi{t), respectively. As noted in the 
previous section, the calculation of these functions requires some care; because the system is not in a statistically steady 
state. Therefore the spectral density of the fluctuations depends on time and the correlation function rxi{t,T) =< 
xi(t—T)xi(t) > depends on both t and on the time difference r. More precisely, the correct definition of the correlation 
function is actually given by 

<Xl{t~ T)xi{t) >-<Xl{t~T)>< Xljt) > 
'^Xi{t-T)0'xi{t) 
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FIG. 7: (a-b) width (in arbitrary units) of the correlation functions of the fluctuation x{t) — xo{t) for the saddle- node model 
(|25|l with a — 10"'^. The correlations functions are defined by equation (|50|l and (|29|) for curves (a) and (b) respectively; (c) 
Spectral width (a. u.) calculated with Wg = to, for a = 10"'^. (d) g{t/tc) for a = 10^'^ (blue), for a = 10~^ (red). 



The latter defiiiition of correlation is not readily accessible in experimental situations, because it requires knowledge 
of the time dependent variance, virhich is difficult to estimate from a single sample. A more accessible tool is often 
used; it is given by the expression 

, , < Xl{t-T)xiit) > - < Xl{t-T) >< Xljt) > 



which coincides with the correct expression if the variance is the same at time t and t — r only. If the variance changes 
noticeably during the time interval of duration r, the latter expression is biased, since T = T' cr(tl^T} ■ expect such 
a discrepancy to manifest itself close to the catastrophe, since the variance increases there by a large amount. To 
illustrate this point we show in Fig|71Ja-b) the correlation functions of the fluctuation x{t) — X(){t) for the solution of 
the saddle-node model, as defined by equations (l50|) and ((29)) respectively. The strong increase of the correlation time 
before the catastrophe is noticeably truncated when using the biased expression (|29p . Using the correct definition of 
the correlation function (|50p , the correlation time increases by a factor of ten over a time interval of order Iq before the 
catastrophe, while the enhancement is only about 3 when using the biased expression (|29p . Moreover the enhancement 
is followed by a drop in the latter case. Nevertheless both curves display well the critical slowing down effect, which 
manifests as a growth of the correlation time close to the transition. The increase of the correlation length before the 
catastrophe is understood by looking at the formal expression ([?7)) . The second derivative of the potential vanishes 
a.t t = to, and that leads to the flatness of g{t) in the time domain < t < tc{a), as shown in FiglTl^d). This time 
domain could therefore be identified as a "precursor time", of order few t^. 

Now let us consider the spectrum of the fluctuations in order to compare with the experimental results. A time 
dependent spectrum can be defined formally by the (real) Wigner transform 

/ — oo 
dr <xi(t-r)xi(t) >e-2*-■^^ (30) 
-oo 

that has to be modified for numerical applications, either by using a filtering procedure like the one used in the 
previous section or by introducing a slipping window. In this section we use the latter process. Choosing a Gaussian 
window function of width Wg , the numerical spectrum is given by 

S,At,^)^{\ rdre-(^)'xi(T)e-2--H, (31) 
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where ta^t are the numerical integration time boundaries. This expression does not take into account the variabihty 
of the variance; therefore we expect the observed change in the spectral properties to be biased, as was the case for 
the correlation function F'. 

The evolution of the spectral width A/ (half- height width) is reported in FiglTl^c) in a range of a few times to around 
the transition. The solution io(t) is hightlightcd by the solid red line. The (c) part of the figure illustrates the same 
effect as the correlation function V , but in the Fourier space. The spectral width decreases noticeably from negative 
time of order few Iq, until the time t ~ 1.5 * ^Oi where it grows, because of the bias due to the variance increase close 
to the burst. The decrease of A/ corresponds to a shift of the spectrum towards low frequencies. The important 
result is that this shift occurs well before the transition. It occurs over a time interval of order few to- This result 
agrees with our experiment where to was estimated to be about 10s and the decrease of the spectral width keeps on 
for about 30s ( see Fig|3] where 10s corresponds to 1300 counts in the abscissa scale). The growth of the fluctuations 
and their shift to lower frequencies can be understood as follows. As the transition approaches, the potential V{x, t) 
becomes flatter and flatter, making weaker and weaker the restoring force toward equilibrium. Therefore, at constant 
noise source, the amplitude of the fluctuations driven by this noise source will grow because the damping is ever 
less efficient. Moreover, the typical time scale for this damping will get increasingly larger because of the decreasing 
stiffness of the potential, thus favoring noise at lower and lower frequencies. 

IV. ANANTHAKRISHNA MODEL 

This section is devoted to an analysis of solutions of a set of equations devised for describing creeping in solids. 
More precisely our purpose is to show that, in a range of parameters this equation exhibits a dynamical saddle node 
bifurcation. As in the "abstract" model of the previous section, this bifurcation is also preceded by a slowing-down 
of the fiuctuations triggered by an external source of noise. Creeping phenomena in real materials are complex and 
difficult to predict quantitatively, despite decades of efforts on theoretical models. We have chosen to consider a model 
developed by Ananthakrishna for creeping in strained solids. To make things simpler, we have only used its version 
without space dependence in the quantities involved. Note that the introduction of space variables would lead to an 
aperiodic creep signal more realistic than the periodic signals of the present model, however we conjecture that it 
should not affect the main result of our study (the emphasis of a precursor signal over a given time interval) . The AK 
model considered here is a set of three coupled non linear ordinary equations with three dimensionless parameters (a, 
h and c, where the letters a and h have no connection with the same symbols used previously). The unknown time 
dependent quantities are three scaled variables corresponding to three density types of dislocations, x{t)^ y{t) and 
z{t), representing respectively mobile, immobile, and those with clouds of solute atoms that mimic Cottrell's idea. 
The model equations write, 

bx{t) = G{x, y) = {l- a)x{t) - hx{tf - x{t)y{t) + y{t). (32) 



y{t) = F{x, y, z) = hx{tf - x{t)y{t) - y{t) + az{t). (33) 



hz{t) = H{z, x) = c{x{t) - z{t)). (34) 

where the variable x{t) stands for the elongation rate. As one can check, if the three variables are positive at 
time zero, they remain so at later times if a, h and c are positive, as assumed. The relative elongation (or strain, 
or creep) L{t) of the solder wire is the time integral of x{t), L{t) = x{t')dt' . Solutions of equations 
have been extensively studied their shape and duration versus the parameter values are given in jl^. Recall 
that relaxation oscillations are depicted for small values of the parameter b only. For c larger than a certain critical 
value Ccr, depending on a, b and not written here, the solution is stable. The fixed point coordinates are given by the 
expressions 

Xs = for a<\ 

a;s = (l-a)(l + V2) for a > i (35) 

at lowest order for the small parameter b. This fixed point becomes unstable by a Poincare-Andronov bifurcation for 
small values of 6, only under the condition that c becomes smaller than c^r ■ In a large range of parameter values 
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the limit cycle associated with the variable x{t) displays relaxation oscillations characterized by slow steps and fast 
bursts. 

In this range of b and c small, we focus on the case c << b, for two reasons. First because it allows the set of 
equations to be reduced to the generic saddle-node equation introduced in the previous section. Secondly because it 
leads to x{t) solutions looking approximately like our experimental data (the strain rate x{t) has to be compared with 
the data files labeled y{n) in sectionjlll). An example of the numerical solution of equations is given in FiglS] 

In this case the burst amplitude drawn in figure (b)) is noticeably larger than the experimental one shown in Fig|2]-(a), 
nevertheless we have chosen the parameter values for pedagogical reasons, in order to give a clear representation of 
the full cycle of the relaxation oscillations. Consider the evolution of the elongation rate x{t). The limit cycle of 
duration T = 180 displays four stages. In the first step, (0 < t < 130, the flow evolves slowly and x{t) takes values of 
order unity. This stage is followed by a fast jump of x{t) at i ^ 130, then by a short slow stage (130 < t < 144) with 
high values of x{t), of order 1/6, finally followed by a fast decrease of x{t) at t ^ 144. 



U, V, teta oraBp(t) 




FIG. 8: (a) Solution of the 3D flow((32l)-(|34l) for a = 0.65 , b = 410"^, and c = 6/100 . The three curves are x{t)/20 (blue), 
y{t) (purple) and z{t)/8 (yellow), (b) Creep signal L{t). 



A. Stability of the 3D flow 



We first consider the linear stability of the 3D flow, especially along the first stage described just above, which 
precedes the fast jump. In the next subsection we derive the local form of the model close to the burst. We shall 
prove that both approaches provide the same information, namely an estimation of the precursor time value. In the 
slow regime, assuming a 3D flow of the form x{t) — xo(t) + SxexpXt (and similar expressions for y(t),z(t)) close to 
the trajectory a:o(^)j yo{t)i zoit), the exponents A are the eigenvalues of the Jacobian matrix 

{l-a-2bxo-yo)/b {l-xo)/b \ 

2bxo-yo ~{xo + l) a , (36) 

c/b -c/b J 

or solutions of the equation 



+ a2A^ + aiA + flo = 0, (37) 

with 02 = {c-{l-a-yo))/b + 3xo + l, ai = {2xo + (ya + a-l) /b){xQ + l + c/b) + {xo + l)c/b- {xo-l)iyo~2bxa)/b, 
and qq — — c(l + a;o)(2&a;o + yo + a — ^) — {yo ^ 2bxo){xo — 1) + a{x — — l)/6^. 

One of the eigenvalues is real and negative all along the trajectory, while the other two become complex conjugates 
in the slow regime, their real part crossing zero before the burst at time tiyap ~ — 3, see Figl9l The time interval 
tc — tiyap depends on the values of the parameters a, &, c. It is generally a small fraction of the limit cycle period. 
Below it is shown that the numerical value of tc — tiyap is nearly equal to the intermediate time scale to for the model. 



B. Normal form close to the burst 



We now consider the behavior of the 3D flow in the vicinity of the burst, and derive the normal form of the AK 
model close to B. 
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FIG. 9: (a) Real and (b) imaginary part of one of the two eigenvalues crossing the real axis before the burst for the case of 
the 3D flow illustrated in figure (1). 

1. slow manifold 

We focus on the first part of the limit cycle, preceding the burst of x(t), and try to understand how the trajectory 
leaves the slow manifold (SM) that we define now. Because c << 6 << 1, the slow stages are described by canceling 
the right-hand side of equations that reduces the 3D flow to the ID flow 



F{x,y,z) = Q 

From the first equation the variables x(t) is an explicit function of y(t), 

x{y) 



(38) 



-{y^l + a) + ^{v-l + aY+Ahy 



Inserting this expression into the second equation 



2h 

allows to also express the variable z{t) in terms of y{t), 



(39) 



z{y) 



^bx{y)^ + y{x{y) + 1) 



(40) 



In the phase space {z,y) expression PO)) defines the slow manifold, which is illustrated in Fig. ([TU)) . red curve. On 
the positive slope parts of the slow manifold the ID flow obeys the differential equation 



y(.t) ^ l — {xiyit)) ^ ziym 

bZ,y 

, which becomes singular at critical points defined (on the SM) by the relation 

Z.y = 0. 



(41) 



The 3D flow (closed blue curve) follows the path A ^ B ^ C ^ D. At the critical point B, the trajectory leaves 
the slow manifold, jumps to point C (fast stage B C), then it follows the portion (C — > I?), and finally returns to 
the SM in A. We consider below the exit of the SM close to the critical point B. 



2. Critical point B 



The coordinates of the critical point B in the phase space x,y,z are solutions of equations ([5^ and (HT|) . From 

dx 
dy 

{x-lf 



equation ([M)) we derive an expression for x^y — || 



1 - a + - 2bx ' 



(42) 
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FIG. 10: SM z{y) (red curve) for a — 0.65 , b = 4.10 , and parametric plot of the 3D flow (blue curve, numerical solution 
z{x{y{t)),y{t)) of equations for c = 6/100. 



which should be identical to the expression of x^y taken from the relation (HJ), 



- 1 



36a;2 - (26 + 1 - a)x' 
These two expressions are identical if xb is a root of the cubic polynomial equation 

x{x - 1)(1 + 26 - 36a;) = (a; + 1)(1 - 26x + bx^). 



(43) 



(44) 



where 6 = 6/(1 — a). The positive solution xb and the corresponding value of yB = ((1 ~ a)xB — bx'^)/{xB ^ 1) and 
zb satisfying equations (j39p -(l40 p respectively are given at leading order with respect to the small parameter 6 by 



xs = x/2 + l + 5±|^T^ 

y5 = (l-a)(l + V2)-^T^ 

^ (3+2v/2)(l-a) _ ,10+7^/2 



(45) 



valid in the parameter range 6 << (1 — a). 

We next derive the normal form describing the dynamical behavior of the solution close to the critical point. 
Assuming that the variable x follows adiabatically the variable y according to equation (j39p . the original system 
reduces to 



y = F{y,z) 

z = iix{y) - z), 



(46) 



defining the dynamics for the portion of the trajectory near B, including the burst B ^ C. This is confirmed by the 
numerics: starting from the point B , with numerical initial conditions (|45p . the behavior of the solutions of equation 
(|46)l agree well with the 3D flow, see Fig. [11] 



)L*Xt,Yt,facz»Zi,Uiic s^ule eq tr 



facx*X-red from B,Y,facz*Z-blue,3eqs from B 




' 5 ' ' ' ' 10 ' ' ' ' 15 ' ' ' ' 20 ' ' ' ' 25 ' ' ' ' 30 ' (b) 

FIG. 11: Solution from point B to point C, (a) Solution of the ID flow equation H34p with relations (|39p -(|40p. (b) solution of 
the 3D flow equations 



Since we are interested in the description of the solution before and at the burst, in the intermediate regime where 
X remains close to a;^ , we shall pursue our analysis by canceling the terms 6a;^ in p9p and in F, because bxB << 1. 
This leads to the system 
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y = F{y,z) (47) 

z = U^iy) - z), 

where F = —xy — y + az. Close to B and during the burst, the variable z is essentially constant, due to the small 
value of the ratio c/b, whereas the variable y jumps toward smaller values, as shown in Fig. 1101 Therefore at leading 
order with respect to the small parameter c/b, the solution of the last equation is given by 

z{t) = ZB + ^{xB ~ ZB)t. (48) 

in the vicinity of B. Inserting this local solution (|48p for z{t) into the system (|T7)) . the exit from the SM is then 
described by a single equation for the local variation 5y — y — yB of y of the variable y , which is of the form 

5y = Y,-AF^y^)B5y'' + l5t, (49) 

n>2 

where 

J = {ac/b){xB - zb), (50) 

and F^"-* sets for the n*'' derivative of F with respect to the variable y, taken at point B ( note that {F^y)B = at 
the critical point B). Using the reduced system (|47|) . for the first two derivatives of F with respect to y we write 

F^y^ — -^yX^y2 — 2x^y ^^^-^ 

Expanding all expressions close to B at leading order with respect to the small parameter 6, one obtains 

{x^y)B = -JZ^ 



_ 4^/2 

y^)B 



(X y2)B = -TT^ 



{Fy^- 
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At this stage we can limit the series expansion in equation (|49p . A rough approximation for the series to converge 
is given by {Fy3)BSy < Z{F^y2)B- Using expressions ([5^ and the first equation (|Tfl) gives the range of variation 



|<5a;| < V2. ^^^^ 
In this range the local form of the AK equations close to B becomes 

5y = \{F^y2)B5y'' +^5t. (54) 

which is identical to equation ([T2|). 

Finally , setting T = {^^^^^^f^fl^St, and Y = [ »^.«^^^^)' ]V^^y, the relation ^ takes the generic form ^ without 
any parameter recently proposed as a possible description of a signal before a catastrophe . 

Y = -Y^ - T. (55) 

As shown in section IIIII the solution of equation displays an intermediate time scale of order 
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which could be used to predict the burst, because the critical slowing down effect occurs during this time interval. In 
the AK equations the intermediate time scale, or precursor time, is given by the relation 



where the parameter 7 is given by ([50)) . Recall that to = tc ~ t^yap was shown to be the time delay between the 
catastrophe and the time t^yap at which the instability builds up along the slow manifold. For the parameter values 
of the above Figures, we have to ^ 3.3, which agrees well with the Lyapunov analysis, see Fig. [51 and also with the 
precursor time deduced from the spectral analysis presented in the next subsection. 

In summary we have proven that the AK model, although formally different from the van der Pol model, has a 
normal form close to the critical point that is consistent with the dynamical saddle-mode model equation studied in 
with an intermediate time scale given by expression ([57p . Consequently, as for the saddle-node model, the AK 
model should display a response to noise with a strong increase of the correlation time occurring with a few times to 
before the burst. 



C. Response to noise before B 

The response to noise of the system ([32|l - p4[) is studied by setting {x — xo{t)+Ub{t), y{t) — yo + Vb{t), z{t) = zo{t)-\- 
Ob{t)), where xo{t),yo{t), zo{t) is the solution of the noiseless AK equations, and the vector V{t) — ui,{t),Vb{t),9b{t) 
characterizes the fluctuations of the response to a noise source. These fluctuations result from the introduction of 
noise terms (either multiplicative or additive) in the original system. In the case of multiplicative noise sources, the 
response of the AK equations is a solution of the system 



xit) = (1/6)[(1 - a)x{t) - bx{tf - x{t)y{t) + y{t)]{l + ej^)) 

m = [bx{tf - x{t)y{t) - y{t) + az{t)]{l + tyfy(t)) (58) 
z{t) - {clb)[x{t)- z{t)]{l + €JS))- 

We report below the result of the numerical study of the correlation function given by expressions (|50p and ([^5[) for 
small amplitude noise. In this case the variance increases slowly along the SM, following approximately the evolution 
of the variable xo{t), as reported in Fig. [T^(a). Figures (b-c) displays the evolution of the correlation function before 
the burst, along the path B ^ C. The width of the correlation function F clearly increases as the burst is approached, 
showing an evident critical slowing down effect. The maximum growth of the correlation time occurs close to the 
time t ^ tiyap where the SM becomes linearly unstable, persisting until the burst. This result is in agreement with 
the experimental result given in section [llj showing a shift of the spectrum towards low frequencies before the burst. 

For noise of additive type the width of the correlation function T{t,T) behaves similarly, as shown in Fig. [T3| As 
pointed out in a previous section, using the expression ([^ to calculate the correlation time of the solution gives a 
biased result close to the burst where the variance is time dependent. In this case the increase of the width before 
the burst is reduced as illustrated in Fig. ITBTb-d). 

In order to present the AK model results by using the same analytical tools as used for the experiment, 
we have calculated the cumulative spectrum of the response Ub(t), Cs{t,i^) = dv' S{t,v') where S{t,v) =< 

I Jj*_^j dt'ub(i') exp^*'^'^*' P > , is the spectrum of the fluctuations for the response signal sampled during the time 
interval {t — At,t), in analogy with equations ©-([31). The characteristic spectral width oj defined by the relation ([5]) 
evolves in time as illustrated in Fig. [Ufor the AK model with additive noise. The horizontal segments corresponding 
to abscissa {t— At, t), have ordinates w, the numerical value of the spectral width calculated during this time interval. 
The decrease of uj close to the burst shows well the expected shift of the spectrum towards low-frequencies, observed 
also in the experiments. 

In conclusion the AK model displays a range of parameter values where one can see the critical slowing-down 
observed also in the experiments. We have studied the Ananthakrishna model for the case of small b and c parameter 
values, with c<^b. In this case the burst occurs close to a critical point, as noticed in [15]. Close to this critical point, 
the 3D flow (namely a set of three coupled ODE's) can be reduced to the 2D system ([47)) with a sort of Langevin-like 
source term (noise term with small amplitude noise)of the form ([^- ([TU)) discussed in section [Till The present system 
differs from the van der Pol equation; however, close to the critical point both models take the form of the dynamical 
saddle-mode model equation studied in 

The Ananthakrishna model displays a response to noise with a strong increase of the correlation time occurring 
close to the burst, or a shift of the spectrum towards low-frequency components. For this model the "precursor time" 
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FIG. 12: Multiplicative noise: (a) variance of the response (red curve) as a function of time, compared to xq (blue curve), 
in Log scale, (b) Correlation functions r(t, r) at times t = 50 -green , t — 100-orange and t = 129-black curve versus r. (c) 
Half-height-half-with of the correlation function T{t, r) versus time t. The input data are those of Figure noise amplitudes 
are ei = 10~^ for the three independent noise sources fi{t). 
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FIG. 13: With additive noise, (a-b) Correlation functions at times t — 50, 100, 129 (dashed, Blue, Red curves respectively), 
(a) r(t, r), (b) T'{t,T) versus r. (c-d) Evolution of the half-width of the correlation functions along the trajectory before the 
burst, (c) half-width of r(t, r) (d) half-width of r'(t, r) both curves drawn versus t. Same parameters and input data as in 
Figures ()12|1 . Noise amplitudes are — 10~*, 10"^, 10^^ for the three independent noise sources fi{t). (e) Evolution of the 
variance along the path (B — )■ C) 




FIG. 14: Cumulative spectral width cj of the response to additive noise Ub{t), versus time before the burst, same data as in 
other figures . 
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is a time interval of order several times tc — tLyap, also equal to several times the intermediate time scale which depends 
on the values of the parameters a, &, c as given by equation A detailed quantitative comparison between the AK 
model and the creep experiment will require a more complete study, which is beyond the scope of the present work. 
In particular it should be noted that including spatial variables in the AK model leads to chaotic solutions, which 
is in better agreement with our experiment where the limit cycle period may indeed vary by a factor of 3 from one 
recorded data set to another. 



V. CONCLUSION 



We have shown that the dynamical model of saddle-node transitions recently proposed to foretell catastrophes is 
applicable to describe the physics of collective dislocations. The experimental signal of the plastic deformation of 
the eutectic mixture of Sn-Pb subjected to a constant stress presented above clearly displays the well-known critical 
slowing down effect, with a precursor time of order 1/10 of the relaxation oscillation period. This observation is 
shown to agree with the response to noise of a physical model proposed by Ananthakrishna for describing creeping in 
ductile materials. This physical model has 3 parameters, therefore it would be quasi-impossible to make a quantitative 
comparison between the experimental data and the theoretical model. Here we show the AK model has a range of 
parameter values for which the slow-fast transition is described by the dynamical saddle-node model mentioned above. 
When adding a small noise to the AK model, we show that it indeed displays the critical slowing down scenario, and 
we derive the expression for the precursor time in terms of the 3 parameter values. This gives a relation between the 
parameter values of the AK model in which the precursor time (i.e. the time during which the most intense spectral 
components clearly shift toward low frequencies) is of order one tenth of the relaxation oscillation period, as observed 
in the experiment. The important result is the following: for such systems which display relaxation oscillations with 
saddle-node transition, we show that there exists an intermediate time scale between the slow and fast regimes, which 
can be used to foretell catastrophes. Unlike the majority of precursor tools that have focused on signals in the time 
domain, we have looked at fluctuations in the frequency domain. These spectral precursors are prone to fewer false 
alarms. 
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